Structures and topological transitions of hydrocarbon films on quasicrystalline surfaces 
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Lubricants can affect quasicrystalline coatings surfaces by modifying commensurability of the 
interfaces. We report results of the first computer simulation studies of physically adsorbed hy- 
drocarbons on a quasicrystalline surface: methane, propane, and benzene on decagonal Al-Ni-Co. 
The grand canonical Monte Carlo method is employed, using novel Embedded Atom Method po- 
tentials generated from ab initio calculations, and standard hydrocarbon interactions. The resulting 
adsorption isotherms and calculated structures show the films' evolution from submonolayer to 
condensation. We discover the presence and absence of the 5- to 6-fold topological transition, for 
benzene and methane, respectively, in agreement with a previsouly formulated phenomenological 
rule based on adsorbate-substrate size mismatch. 



Friction can become vanishingly small between incom- 
mensurate interfaces (superlubricity) [l|, Q ■ The aperi- 
odic structure of quasicrystals surfaces makes them the 
ideal candidates for this phenomenon. Indeed, the special 
frictional properties of quasicrystals have been confirmed 
in a series of experiments using atomic force microscopy 
on single grain quasicrystals [1, 0, H, H, H HI • Although 
its origin is not yet completely understood, the evidence 
suggests that poor coupling of phonons at the interfaces 
may play a major role @, |3i S]- In addition, there is 
also evidence that oxidation decreases the friction even 
further 0,0]. In fact, even before these single crystal ex- 
periments were performed, experiments on quasicrystal 
coatings in air did show low- friction behavior [ld|, lead- 
ing to applications involving moving machine parts and 
non-stick cookware [llj. 

Since the low friction of quasicrystals is clearly re- 
lated to their structure, the interaction of the lubricant 
with the quasicrystal and the structure of the lubricat- 
ing film is particularly important. In most applications 
involving machine parts, additional lubricant would be 
needed to address the frictional contributions of grain 
boundaries and asperities. Studies of metal and rare 
gas adsorption on quasicrystals indicate that both pe- 
riodic and aperiodic structures can occur in thin films 
13, E mill, El, [l3- However, little is known about 
the interaction of hydrocarbons with quasicrystal sur- 
faces 

In the present manuscript, we explore the effect of 
structural and symmetry mismatches on the ordering of 
hydrocarbons by evaluating the nature of hydrocarbon 
adsorption (methane, propane, and benzene) on a qua- 
sicrystalline decagonal surface, namely, the tenfold sur- 
face of AlysNiioCoiv. The simulations are performed us- 
ing grand canonical Monte Carlo method (GCMC), with 
which we have extensive experience on smooth [20l . I21I [2^ 
and corrugated substrates H, 17, H, 0,1231. Using 



the GCMC method, we compute the adsorption proper- 
ties for specified thermodynamic conditions. We take a 
tetragonal unit cell, of side 5.12 nm, with a hard wall at 
10 nm above the surface to confine the coexisting vapor 
phase without causing capillary condensation (sufRcient 
to contain 25 layers of benzene). We assume periodic 
boundary conditions which, although sacrificing the ac- 
curacy of the long range QC structure, do not interfere 
with short-to-moderate length scales, representative of 
the hydrocarbon order. The substrate is reproduced with 
an 8-layer Al-Ni-Co slab, where the atom coordinates 
are derived from an experimental low-energy electron 
diffraction (LEED) study [26]. The inter-molecular in- 
teractions (adsorbate-adsorbate) are calculated as a sum 
of pair interactions between atoms. Buckingham poten- 
tials are used for methane [27,, 28] and benzene 2^ 3^ 
while a Morse potential is employed for propane |3l| . 
We have developed embedded-atom method (EAM) po- 
tentials (3^ to model the many-body intra-molecular in- 
teractions, adsorbate-substrate (C-Al, C-Co, C-Ni, H-Al, 
H-Co, H-Ni), and substrate-substrate (Al-Al, Al-Co, Al- 
Ni, Co-Co, Co-Ni, Ni-Ni). The EAM embedding func- 
tions are taken as natural cubic splines using charge den- 
sity functionals from Herman [33|, while the EAM pair 
energies have Morse potential forms [34| with Haftel's 
mixing scheme [s^, if necessary. The EAM potential 
parameters are fitted from ab initio energies calculated 
using the Vienna Ab-initio Simulation Package (VASP) 
(36j with exchange-correlation functionals as parameter- 
ized by Perdew, Burke, and Ernzerhof (PBE) [37[ for the 
generalized gradient approximation(GGA), and projec- 
tor augmented-wave (PAW) [11] pseudopotentials. The 
EAM potential parameters, fitted using the SIMPLEX 
method ^], are summarized in the EPAPS material [ioj. 
Quasicrystal approximants (crystals having similar short- 
range order) are used to address the periodic bound- 
ary conditions required by the ab initio package: we use 
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Al29Co4Ni8 (dBl), AlivCosNia (dHl), and Al34Co4Nii2 
(dH2) [4l[. Our calculations show that alkanes and ben- 
zene do not dissociate on such substrates, which do not 
undergo any considerable relaxation upon the adsorption 
of the molecules. Thus, as a first approximation, the 
substrate and the molecules can be considered rigid, al- 
though the molecules are allowed to explore all rotational 
degrees of freedom (42j . 




FIG. 1: (color online). Computed potential energies for 
(a) methane and (b) benzene on Al-Ni-Co, obtained by 
minimizing V{x, y, z, 9, 4>, ip) of a single molecule with re- 
spect to {z,^,9,<j)), variations. The average values are 221 
meV/methane and 931 meV/benzene. 

The adsorption potential of a single molecule, 
Vminix,y), calculated as the maximum depth as a func- 
tion of normal coordinate (z) and Euler angles (6', 0, -0), 
is both deep and extremely corrugated. Figure [T] shows 
Vminix,y) for methane and benzene: the dark spots in- 
dicate strong binding sites. The average adsorption en- 
ergy is 221 meV/methane, 374 meV/propane, and 931 
meV/benzene. 

The symmetry of the adsorption potentials for 
methane and benzene reflect the pentagonal symmetry 
of the substrate, as illustrated in Fig. [T] Propane follows 
a similar trend as methane, with somewhat less corruga- 
tion due to its larger size. The most attractive adsorp- 
tion positions for methane and propane [40| are located 
at the centers of pentagonal hollows having five Al atoms 
at the vertices. Conversely, the most attractive sites for 
benzene are the Al-centered pentagons with 3 Al and 2 
Ni atoms at the vertices. These hollow and Al-centered 
pentagons alternate every 36° around the z axis. 

By simulating the adsorption of a single molecule, a 
general trend is observed for the orientation of the ad- 
sorbant: the smaller the molecule, the more variation in 
the rotation of the ground state. Methane's ground state 
is highly degenerate, propane adsorbs with its axis form- 
ing a small angle with respect to the substrate (the angle 
varies from 0° to 10° depending on the adsorption site). 
Benzene adsorbs with its plane parallel to the substrate. 

Figure [2] shows the computed adsorption isotherms 




FIG. 2: (color online). Isothermal adsorption densities (pjv 
in molecules/nm'^) of hydrocarbons on a decagonal Al-Ni-Co: 
(a) methane, (b) propane, and (c) benzene. The simulation 
temperatures are reported in blue (in Kelvin) . The insets rep- 
resent the densities along the z direction at P corresponding 
to points "d" . 

Pn (densities in molecules /nm^) at different tempera- 
tures for methane (T=68, 85, 136, 185 K), propane (T = 
80,127,245,365 K), and benzene (T = 209, 270, 418, 555 
K) as functions of pressure P. The plotted quantities 
are the thermodynamic excess coverages (differences be- 
tween the total number of molecules and the number that 
would be present if the cell were filled with uniform va- 
por). The substrate is very attractive (Fig. [T]) hence 
the hydrocarbons experience complete wetting up to the 
highest temperature close to the critical temperature. 
At low temperatures, the formation of the first layers 
is evident from the first step in each plot (the most left 
isotherms) . The formation of further layers is not clearly 
observed in contrast with the observations for noble gases 



15, [l^, 17, 24 1 . Nevertheless, layering in the adsorbed 
film is revealed by the insets of Fig. [2] showing the ad- 
sorption densities along the z direction at the pressures 
corresponding to points "d" of the lowest temperature 
isotherms. The relative positions of zi, Z2, and Z3 indi- 
cate that methane adsorbs mostly with three hydrogens 
anchored to the substrate, while propane adsorbs mostly 
with five hydrogens near the substrate (two from each 
end and one from the middle segment). 

In the submonolayer regime (point "a" in Figure [2^ a)), 
methane adsorbs preferentially at the strong binding 
sites. By increasing P, a methane monolayer forms with 
pentagonal ordering commensurate with the substrate 
(point "c" in Figure a)). Similar configurations are ob- 
served at all studied temperatures. An example of such 
ordering is illustrated in Fig. [3l^a) at T=68K: the den- 
sity profile p(x^ y) and the Fourier transform FTcm of the 
density of the center of mass confirm pentagonal order- 
ing (10 discrete spots representing 5-fold axes in FTcm). 
Similar plots for propane at 80 K and benzene at 209 K 
are shown in Figure [3] (b) and (c) , respectively. At all 
simulated temperatures, benzene has 6-fold order as in- 
dicated by its FTcm characteristic of triangular lattice. 
Unlike methane and benzene, which adsorb in a well- 
defined structure, propane forms a poorly-ordered 5-fold 
arrangement (clusters of molecules which form pentagons 
can be seen in the EPAPS material ji^l). The FTcm in- 
dicates a distortion involving a compression along one 
direction. This is due to the interplay between the linear 
form of propane molecule and the quasicrystal structure. 
Figure |4] shows a histogram of the orientations of the 
propane's axis projected on the xy-plane (panel (a)) and 
the a:z-plane (panel (b)), corresponding to the density 
plot in Figure [Sjb). Indeed, propane adsorbs with its 
axis parallel with the substrate (up to 10°) and preferen- 
tially oriented along y-axis. 

Figure m^a) illustrates a superposition of the methane 
monolayer at 85 K and the top layer of substrate atoms. 
Methane molecules are present in every hollow pentagon 
defined by five Al atoms. These pentagons, correspond- 
ing to the dominant binding sites and depicted as dark 
spots in Fig. [l]^a) , are responsible of stabilizing the qua- 
sicrystalline structure of the methane monolayer. A sim- 
ilar plot for benzene at 209 K is given in Figure [5Kb). A 
domain having 6-fold order can be observed in the left 
part of the figure. Dotted circles corresponding to the 
benzene molecules, along with two red lines illustrate the 
orientation of this domain, which makes an angle of 18° 
from the y-axis. Thus, the benzene lattice is oriented 
along a symmetry direction of the quasicrystal, as ob- 
served previously for the Xe monolayer. 

Previosuly, by studying adsorbed monolayer of rare 
gases 16|, |l7| , we found that the crucial parameter in de- 
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FIG. 3: Monolayer density profiles p(x, y) and Fourier trans- 
forms of the density of the center of mass FTcm for: (a) 
methane at 68 K (point "c" in Figure [2^ a)), (b) propane at 
80 K (point "c" in Figure [SKb)), and (c) benzene at 209 K 
(point "c" in Figure [JJc)) adsorbed on decagonal Al-Ni-Co. 



(a) 

E 



(b) 

E 



orientation on XY plane 


< along y-axis 

J aiu .aniui.. ,.iu ini i ah 


I orientation on XZ plane 

II <— horizontal horizontal 1 

ML 1 1 HIT 



15 30 45 



75 90 105 120 135 150 165 180 



FIG. 4: Histogram of orientation of propane's axis on (a) 
XY plane and (b) XZ plane at monolayer coverage at 80 K 
adsorbed on decagonal Al-Ni-Co. The corresponding density 
plot is depicted in Figure (Sjb). 

termining the overlayer structure is the relative size mis- 
match between adsorbate and substrate's characteristic 
length. The mismatch is defined S'S (5m = (^r — '^r)/'^r 
where dr is the distance between rows in a 2D close- 
packed arrangement of adsorbates 



17, m 



and Ar 

0.381 nm is the quasicrystal's characteristic length [Ij 
Thus, 5m measures the relative mismatch between an ad- 
sorbed {111} closed packed plane of adsorbates and the 
QC surface. The proposed rule states that the ordering 
transition exists if and only if (5m > 17[ . For methane, 



by starting from the lattice parameter of the cubic cell 
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M l45|, we have SmiCUs) ^ -0.055. For benzene, by 
using dr obtained from the simulation of the 6-fold phase 
(using a flat substrate with the well depth of 931 meV), 
we get SmiCQHe) ~ 0.617. Both methane and benzene 
satisfy the mismatch rule extending its applicability (the 
transition can not be defined for propane since it does 
not form pentagonal or triangular arrangements). 



Teragrid Partnership (Texas Advanced Computing Cen- 
ter, TACC) for computational support (MCA-07S005). 
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FIG. 5: (color online), (a) methane monolayer at 85 K (b) 
benzene monolayer at 209 K adsorbed on decagonal AI-Ni-Co. 
Top layer of substrate atoms are plotted in blue (Al), green 
(Ni), and red (Co). 

TABLE I: Summary of adsorbed rare gases and selected hy- 
drocarbons undergoing (or not) the commensurate incom- 
mensurate transition on-Al-Ni-Co. iNe, iXe, and dXe are 
inflated/deflated noble gases [l^]. Triangles and pentagons 
indicate triangular lattice and fivefold structure, respectively. 





S,n 


transition monolayer 


methane 


-0.055 


No 


O 


benzene 


0.617 


Yes 


A 


Xe [15, 


0.016 


Yes 


A 


Ne [IT] 


-0.311 


No 


A + O 


Ar [IT] 


-0.158 


No 


A + O 


Kr [ir| 


-0.108 


No 


A+ O 


iNe^i' [IT] 


0.016 


Yes 


A 


dXe<i',dXe(2) [17] 


-0.311,-0.034 


No 


A + O 


iXe(i),iXe(2) [17] 


0.363, 0.672 


Yes 


A 



Table U summarizes the results obtained for all hydro- 
carbons and rare gases that we have studied so far. When 
the transition is not present, the monolayer does not form 
a structure having long-range periodic order. From the 
point of view of the potential for superlubricity, all of 
the gases that are smaller than Xe do not form periodic 
films. In addition, the linear molecule propane, which is 
larger than Xe, does not. This suggests that the shape 
as well as the size is an important factor in the ordering, 
and it may bode well for the use of linear hydrocarbons 
as lubricants on quasicyrstalline films. 
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